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ABSTRACT 

Intracluster stellar populations are a natural result of tidal interactions in galaxy clusters. Measuring 
these populations is difficult, but important for understanding the assembly of the most massive 
galaxies. The Coma cluster of galaxies is one of the nearest truly massive galaxy clusters, and is host 
to a correspondingly large system of globular clusters (GCs). We use imaging from the HST/ACS 
Coma Cluster Survey to present the first definitive detection of a large population of intracluster GCs 
(IGCs) that fills the Coma cluster core and is not associated with individual galaxies. The GC surface 
density profile around the central massive elliptical galaxy, NGC 4874, is dominated at large radii by 
a population of IGCs that extend to the limit of our data (R < 520 kpc). We estimate that there are 
47000 ± 1600 (random) l^oo (systematic) IGCs out to this radius, and that they make up ~ 70% of 
the central GC system, making this the largest GC system in the nearby Universe. Even including 
the GC systems of other cluster galaxies, the IGCs still make up ~ 30-45% of the GCs in the cluster 
core. Observational limits from previous studies of the intracluster light (ICL) suggest that the IGC 
population has a high specific frequency. If the IGC population has a specific frequency similar to 
high-SV dwarf galaxies, then the ICL has a mean surface brightness of \iv ps 27 mag arcsec -2 and 
a total stellar mass of roughly 10 12 A4© within the cluster core. The ICL makes up approximately 
half of the stellar luminosity and one-third of the stellar mass of the central (NGC4874+ICL) system. 
The color distribution of the IGC population is bimodal, with blue, metal-poor GCs outnumbering 
red, metal-rich GCs by a ratio of 4:1. The inner GCs associated with NGC 4874 also have a bimodal 
distribution in color, but with a redder metal-poor population. The fraction of red IGCs (20%), and 
the red color of those GCs, implies that IGCs can originate from the halos of relatively massive, L* 
galaxies, and not solely from the disruption of dwarf galaxies. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: clusters: individual: Coma — 
galaxies: halos — galaxies: evolution — galaxies: star clusters: general - globular 
clusters: general 



1 Based on observations with the NASA/ESA Hubble Space Tele- 
scope obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astron- 
omy, Inc., under NASA contract NAS 5-26555. 

2 Department of Astronomy, Peking University, Beijing 100871, 
China; peng@bac.pku.edu.cn 

3 Kavli Institute for Astronomy and Astrophysics, Peking Uni- 
versity, Beijing 100871, China 

4 Space Telescope Science Institute, 3700 San Martin Drive, Bal- 
timore, MD 21228, USA 

5 Department of Physics and Astronomy, Johns Hopkins Univer- 
sity 3400 N. Charles St., Baltimore, MD 21228, USA 

6 Department of Physics, University of Durham, South Road, 
Durham DH1 3LE, UK 

7 Dept. of Physics & Astronomy, San Francisco State University, 
1600 Holloway Avenue, San Francisco, CA 94132, USA 

8 National Research Council of Canada, Herzberg Institute of 
Astrophysics, 5071 West Saanich Road, Victoria, BC V9E 2E7, 
Canada 

9 Departmento de Astronomi'a y Astrofi'sica, Pontificia Universi- 
dad Catolica de Chile, Av. Vicuna Mackenna 4860, 7820436 Macul, 
Santiago, Chile 

10 Astrophysics Research Institute, Liverpool John Moores Uni- 
versity, Twelve Quays House, Egerton Wharf, Birkenhead CH41 
1LD, UK 

11 Instituto de Astrofisica de Canarias, 38200 La Laguna, Tener- 
ife, Spain 

12 Department of Physics, Engineering Physics and Astronomy, 
Queen's University, Kingston, ON K7L 3N6, Canada 



1. INTRODUCTION 

1.1. Intracluster Stellar Populations and Hierarchical 
Galaxy Formation 

Massive elliptical galaxies at the centers of galaxy 
clusters — often brightest cluster galaxies (BCGs) and 
sometimes cD galaxies — generally have little ongoing star 
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formation with only minor evolution since z ~ 1, and 
only a shallow relationship between their stellar mass and 
their host cluster mass (e.g., Lin & Mohr 2004; Whiley 
et al. 2008). In the standard hierarchical paradigm, how- 
ever, the most massive halos should be the last to assem- 
ble, and so these galaxies have traditionally presented 
problems for formation models. 

Recent simulations suggest that this paradox can be 
resolved in a picture where the stars that end up in the 
most massive galaxies form early, energy feedback from 
supernovae and active galactic nuclei subsequently sup- 
press star formation, and the assembly of these galaxies 
through dry mergers continues right up to the present 
day (De Lucia & Blaizot 2007; although sec Bildfell et al. 
2008 for evidence that feedback is not 100% efficient). 
Thus, even if there is little star formation at late times, 
these galaxies are still expected to further assemble stel- 
lar mass at z < 2. This predicted increase in the masses 
of BCGs over time, however, may be in conflict with ob- 
servations that show little mass evolution of BCGs from 
z ~ 1.5 to the present (Collins et al. 2009). 

It is expected that the dry merging or tidal stripping 
of satellites should not only contribute stars to the cen- 
tral galaxy itself, but also to an intracluster component 
that has previously been associated with the extended 
stellar envelopes of cD galaxies (Matthews, Morgan & 
Schmidt 1964), and is sometimes labeled as a "diffuse 
stellar component" (Monaco et al. 2006) or simply "intra- 
cluster light" (ICL). This component can make up a large 
fraction of the total luminosity at the center of galaxy 
clusters (Ocmler 1976), and if added to the stellar mass 
of central cluster galaxies, might naturally explain cur- 
rent contradictions between simulation and observation. 
Purcell, Bullock, & Zentner (2007) simulated the forma- 
tion of the ICL from the shredding of satellite galaxies, 
finding that in massive clusters, the ICL can dominate 
the total stellar mass of the combined ICL+BCG sys- 
tem, which is consistent with observations of low redshift 
clusters (Gonzalez, Zabludoff, & Zaritsky 2005; Seigar, 
Graham, & Jerjen 2007). 

In fact, there is increasing observational evidence 
that a significant fraction (10-40%) of the total stellar 
light in a galaxy cluster is intcrgalactic. Starting with 
Zwicky (1951), many detections of low surface brightness 
starlight in galaxy clusters — both in cD envelopes and 
in the regions between galaxies — support the existence 
of substantial intracluster stellar populations (Welch & 
Sastry 1971; Uson et al. 1991; Vilchez-Gomez et al. 1994; 
Gregg & West 1998; Trentham & Mobasher 1998; Feld- 
meier et al. 2002, 2004a; Lin & Mohr 2004; Adami et al. 
2005; Zibetti et al. 2005; Mihos et al. 2005; Gonzalez 
et al. 2005; Seigar et al. 2007; Gonzalez, Zaritsky, & 
Zabludoff 2007; Krick & Bernstein 2007). In nearby 
clusters, there have also been direct detections of in- 
tergalactic red giant branch stars (Ferguson, Tanvir, & 
von Hippel 1998), asymptotic giant branch stars (Dur- 
rell et al. 2002) planetary nebulae (Theuns & Warren 
1997; Mendez et al. 1997; Feldmeier, Ciardullo, & Ja- 
coby 1998; Feldmeier et al. 2004b; Okamura et al. 2002; 
Arnaboldi et al. 2004; Gerhard et al. 2007; Arnaboldi 
et al. 2007, Castro-Rodriguez et al. 2009; Doherty et al. 
2009), novae (Neill, Shara, & Oegerle 2005), and super- 
novae (Gal- Yam et al. 2003). Detections of intergalactic 
light have also been made in compact groups (Da Rocha 



& Mendes de Oliveira 2005). These studies tend to show 
that the richer and more massive the group or cluster, 
the larger the fraction of intergalactic light. 

1.2. Intergalactic Globular Clusters 

Another important clue to the formation of massive 
ellipticals is that those residing at the centers of galaxy 
clusters often host extremely large populations of glob- 
ular clusters (GCs). The star forming events that form 
globular clusters will mostly form stars that end up in the 
field, so it is natural that the number of GCs in a galaxy 
should roughly scale with that galaxy's stellar luminosity 
or mass. However, the ratio of GCs to starlight — usually 
characterized as the specific frequency, Sn (Harris & van 
den Bergh 1981) — has long been known to vary across 
galaxy mass and morphology, with giant central ellipti- 
cal galaxies harboring the largest GC systems and having 
some of the highest specific frequencies. This abundance 
of GCs in galaxy clusters appears explainable if the num- 
ber of GCs scales with either total baryonic mass at the 
cluster center, including hot gas (McLaughlin 1999), or 
the total dynamical mass of the cluster (Blakeslee et al. 
1997). Blakeslee (1997, 1999) observed that the number 
of GCs in galaxy clusters was directly related to cluster 
mass, but the relatively constant BCG luminosity thus 
led to high 5jv- It is possible that the high specific fre- 
quencies are because measurements of the galaxy lumi- 
nosity typically do not include a substantial ICL com- 
ponent, and that the high Sn in central cluster galaxies 
would be more normal if the ICL was included. Galaxy 
specific frequency also varies with galaxy stellar mass (or 
luminosity) in a way that is consistent with the expected 
variation in galaxy stellar mass fraction (or mass-to- light 
ratio) (Peng et al. 2008; Spitler et al. 2008). 

This connection between GCs and total mass has inter- 
esting implications, particularly in massive galaxy clus- 
ters where the predicted build up of stellar mass in cen- 
tral galaxies should be paralleled by the build up of a 
large GC system. If much of the stellar mass in galaxy 
clusters resides in the low surface brightness ICL then 
there should also be a corresponding population of intr- 
acluster GCs (IGCs) that are not gravitationally bound 
to individual galaxies, but directly to the cluster itself. 
Moreover, the detection of point source IGCs in the near- 
est clusters is a much easier observational endeavor than 
measuring the faint ICL, giving us a window onto the 
nature of the diffuse stellar content. 

There are other reasons to expect substantial popu- 
lations of IGCs. West et al. (1993) proposed that GC 
formation may be biased toward the largest mass over- 
densities, i.e. galaxy clusters. West et al. (1995) also pro- 
posed that populations of IGCs were responsible for the 
high Sn seen in cD galaxies. More recently, spectroscopy 
of ultra-compact dwarfs (UCDs) and massive GCs (also 
dubbed dwarf- globular transition objects; Hasegan et al. 
2005) have uncovered a population of compact stellar sys- 
tems in galaxy clusters resembling the most massive GCs 
or dE nuclei stripped of their host galaxies (Drinkwater 
et al. 2003; Hilker et al. 2007; Mieske et al. 2008; Gregg 
et al. 2009; Madrid et al. 2010; Chiboucas et al. 2010). 
These objects, while generally more massive than typical 
GCs and consequently may have different origins, might 
be the so-called "tip of the iceberg" for a large population 
of free-floating, normal globular clusters. 
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In fact, a number of extragalactic GC studies over the 
past few years have strongly suggested the presence of 
IGCs in nearby galaxy clusters. In the Virgo and Fornax 
Clusters, serendipitous discoveries of GCs in intergalac- 
tic regions using HST imaging (Williams et al. 2007), 
ground-based imaging (Bassino et al. 2003) and spec- 
troscopy (Bergond et al. 2007) point to the existence of 
IGC populations. However, it is often unclear whether 
these GCs are truly intergalactic, or are part of the ex- 
tended halos of cluster galaxies (c.f. Schuberth et al. 
2008). A recent study by Lee, Park & Hwang (2010), 
however, used data from the Sloan Digital Sky Survey 
(SDSS) and found statistically significant detections of 
GC candidates throughout the Virgo cluster. In more 
distant galaxy clusters, candidate IGC populations have 
been identified as point source excesses in HST imaging 
(Jordan et al. 2003; West et al. 2011). 

It is possible that some IGCs and intracluster stars 
formed in situ, i.e. in cold, intergalactic gas that never 
accreted onto or was stripped from galaxies. It is also 
possible that IGCs formed very early and at high effi- 
ciencies in dwarf-sized subhalos (e.g., Moore et al. 2006; 
Peng et al. 2008), and whose host galaxies were subse- 
quently tidally destroyed by interactions with the cluster 
potential. Another possibility is that IGCs were formed 
in larger galaxies and were stripped through tidal inter- 
actions with the cluster potential or with other galaxies 
(see, e.g., simulations of Yahagi & Bekki 2005; Bckki & 
Yahagi 2006). 

The formation of the IGC population is obviously 
linked to that of the ICL, although the observed proper- 
ties of the two populations may be different. For exam- 
ple, the detectability of the ICL is highly dependent on 
its surface brightness, whereas IGCs are detectable even 
in isolation. Simulations by Rudick et al. (2009) show 
that the ICL is supplied by tidal streams that originally 
have relatively high surface brightness but then disperse 
to become fainter and harder to detect. ICL studies using 
surface photometry are thus more sensitive to recent dis- 
ruptions, whereas the IGC population is a less temporally 
biased tracer of the full intracluster stellar population. 

The study of extragalactic GC systems has been trans- 
formed by the high spatial resolution imaging of the 
Hubble Space Telescope (HST), with observations of hun- 
dreds of GC systems now in the archives and published 
literature (e.g., Jordan et al. 2009). HSTs deep sen- 
sitivity to compact or unresolved sources, and its abil- 
ity to distinguish background galaxies from likely GCs, 
makes it an ideal tool for extragalactic GC studies. How- 
ever, the relatively narrow field of view of HSTs cam- 
eras and the close proximity of the galaxies being stud- 
ied (D < 100 Mpc) means that most HST studies have 
focused on GC systems directly associated with galaxies. 
Observations of wider fields are usually conducted with 
ground-based telescopes (e.g. McLaughlin 1999; Bassino 
et al. 2006; Rhode et al. 2007), that gain area at the cost 
of spatial resolution. 

1.3. The Coma Cluster of Galaxies 

The Coma Cluster of galaxies (Abell 1656) is one of the 
nearest rich, dense clusters, and is a fundamental target 
for extragalactic studies. Studies of GC systems in Coma 
from the ground using surface brightness fluctuations 
(Blakeslee et al. 1997; Blakeslee 1999) and using HST 



(Kavelaars et al. 2000; Harris et al. 2000; Harris et al. 
2009) all point to large GC systems around the cluster's 
giant elliptical galaxies, particularly around the central 
massive galaxy, NGC 4874 (Harris et al. 2009). Although 
many photometric studies support the existence of an in- 
tracluster stellar light component in Coma (e.g., Zwicky 
1951; de Vaucouleurs & de Vaucouleurs 1970; Welch & 
Sastry 1971; Kormendy & Bahcall 1974; Matilla 1977; 
Mclnick, White & Hoessel 1977; Thuan & Kormendy 
1977; Gregg & West 1998; Calcaneo-Roldan et al. 2000; 
Adami et al. 2005), and even velocities for intracluster 
PNe have been measured (Gerhard et al. 2007; Arnaboldi 
et al. 2007), evidence for or against the existence of IGCs 
is much more muddled. A search for IGCs in Coma using 
ground-based data by Marm-Franch et al. (2002, 2003) 
did not find a surface brightness fluctuation signal that 
would have hinted at the presence of IGCs, mainly be- 
cause of the shallow limiting magnitude of their photom- 
etry. 

At a distance of 100 Mpc, the value we adopt for this 
paper (m - M = 35, Carter et al. 2008), 1' on the sky 
subtends 29 kpc in the cluster, and the mean of the GC 
luminosity function (GCLF) in giant ellipticals is Iyega — 
26.44 mag. This regime of projected areal coverage and 
GC apparent brightness makes it reasonable to conduct 
a contiguous survey of the Coma cluster core for GCs, 
unbiased by the locations of individual galaxies, using the 
HST Advanced Camera for Surveys (ACS) Wide Field 
Channel (WFC). 

The HST/ A CS Coma Cluster Survey is a Treasury Sur- 
vey originally approved for 164 orbits. One of the main 
components of this survey was a contiguous ACS/ WFC 
mosaic of the core of the Coma cluster, making it the 
ideal data set to investigate the existence of intergalactic 
GCs. Hints of this population have already appeared in 
studies of UCDs by Madrid et al. (2010) and Chiboucas 
et al. (2010). This paper presents the first compilation 
and description of IGCs in the Coma cluster. 

2. OBSERVATIONS AND DATA 

2.1. Imaging Data 

The data used in this study are from the HST/ ACS 
Coma Cluster Survey. The survey observations and data 
reduction are described in detail by Carter et al. (2008), 
the catalog generation for the public data release is de- 
scribed by Hammer et al. (2010), and an in depth anal- 
ysis of galaxy structural parameters and completeness is 
presented in Balcclls et al. (2010). We summarize the 
relevant information here. 

The Coma Cluster Survey, as originally designed, con- 
sisted of a large central ACS mosaic of the Coma cluster 
core, and 40 targeted observations in the outer regions 
of the cluster. The central mosaic was designed to be 42 
contiguous ACS/WFC pointings in a 7 x 6 tiling config- 
uration, and covering an 21' x 18' area. Each pointing 
is observed in two filters, F475W (g) and F814W (/), 
with exposure times of 2560s and 1400s, respectively. 
Unfortunately, the failure of the ACS/WFC (January, 
2007) meant that only 28% of the survey was completed: 
19 pointings in or around the central mosaic, and 6 in 
the outer regions. Within the core, the central galaxy 
NGC 4874 was observed, but the other giant elliptical, 
NGC 4889, was not imaged before the ACS failure. De- 
spite the shortfall, the current observations still provide 
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Fig. 1. — (a) The ACS F814W image of the central pointing (Visit 19) containing NGC 4874 and many other bright elliptical galaxies. 
North is left, and East is down. The image is 202" (98 kpc) on a side, (b) The same pointing but after our iterative galaxy subtraction. 
Residuals at galaxy centers are still visible at this contrast level, but the overall large scale gradients in the background light have been 
removed. 



the largest set of deep, high resolution imaging avail- 
able for this important galaxy cluster. Recent studies of 
compact galaxies in Coma (Price et al. 2009) and spec- 
troscopy of Coma cluster members (Smith et al. 2009) are 
part of a concerted effort to study galaxy evolution in the 
Coma cluster built around this HST Treasury survey. 

The ACS data reduction was performed using a dedi- 
cated Pyraf/STSDAS pipeline that registered and com- 
bined images while performing cosmic-ray rejection. 
The dithered images were combined using Multidrizzle 
(Koekemoer et al. 2003), which uses the Drizzle algo- 
rithm (Fruchter & Hook 2002). For this study, we used 
the data drizzled for the ACS Coma Survey Data Release 
2 (DR2). However, except for Visits 3, 10, and 57, the 
F814W images on which the bulk of this paper is based 
are identical to those in Data Release 1 (DR1). 

2.2. Object Catalogs and Galaxy Subtraction 

Our images of the Coma cluster reveal a striking 
amount of detail: cluster members across the mass spec- 
trum, globular clusters, background galaxies, and a few 
foreground stars. Given the different spatial scales of 
these objects on the sky, it is important to generate cat- 
alogs with detection parameters optimized for the sub- 
ject under study. For our purpose of studying the glob- 
ular clusters between galaxies, we used the well-tested 
Source Extractor software package (Bertin & Arnouts 
1996) with parameters optimized for point source detec- 
tion, and which are effectively identical to those used by 
Hammer et al. (2010) in the public data release 1 . Pho- 
tometry was put on the AB magnitude system using the 
zeropoints of Siranni et al. (2005). All magnitudes in this 
paper are AB unless otherwise specified. 

For most visits, the area between galaxies is much 
larger than that occupied by galaxies and the catalogs 
can be considered effectively complete to the same level 

For details, visit the Coma Cl uster Survey website at 
http : //astronomy . swin. edu. au/coma/ 



except in the close vicinity of cluster galaxies. The one 
exception is Visit 19, which contains NGC 4874 and 
many other ellipticals in the cluster core (Figure [TJl) . 
This pointing is nearly entirely dominated by the light 
from one galaxy or another; it is also the one with 
the highest concentration of GCs. To address this, we 
implemented an iterative galaxy subtraction algorithm 
that produced a fully background subtracted image (Fig- 
ure Q]d). 

We subtracted the 10 brightest galaxies from Visit 
19. We started from the brightest (NGC 4874 itself) 
and worked to the faintest of the ten. In each case, 
we first manually masked all the bright galaxies except 
for the one being subtracted. We then used the IRAF 
ellipse and bmodel tasks (Jedrzejewski 1987) to model 
the isophotes of the object galaxy. Because the ellipse 
fitting only occurs out to a finite radius, the resulting 
model will have finite extent, and the subsequent sub- 
traction will leave a sharp discontinuity in the image. 
For convenience of object detection, we extended the 
ellipse-generated models with a power-law fit to the 
last five data points in the profile at fixed ellipticity. 
This allowed for a smooth subtraction out to the bor- 
ders of the ACS image. After subtracting one galaxy, we 
then repeated the process with the next brightest galaxy 
on the subtracted image. After the last galaxy was sub- 
tracted, we used Source Extractor to create and subtract 
a background map that removed large scale variations. 
This last step is important because it allows us to recover 
from any large scale over- or under-subtractions due to 
mismatches between the power-law extensions and the 
true surface brightness profiles of the galaxy. A similar 
technique was used with success by Jordan et al. (2004) 
in ACS images of Virgo cluster galaxies, although that 
was only for single galaxies. 

After galaxy subtraction, we generated catalogs with 
Source Extractor, using variance maps that accounted 
for the extra Poisson noise expected from the subtracted 
galaxy light. These catalogs contain objects much closer 
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to the centers of galaxies, and to NGC 4874 in partic- 
ular. Photometry was obtained by using 3 pixel radius 
(0'.'15) circular apertures with aperture corrections and 
zeropoints from Sirianni et al. (2005). Unless otherwise 
specified, all magnitudes in this paper are on the AB sys- 
tem. These objects are included in the DR2 catalogs of 
Hammer et al. (2010). 

2.3. Completeness 

We use artificial star tests to quantify the spatially 
varying detection efficiency across our images. This is 
particularly important for the galaxy-subtracted image 
containing NGC 4874, where the bright galaxy light af- 
fects the depth of our observations. 

We first use routines in DAOPHOT II (Stetson 1987) 
to construct an empirical PSF using bright point sources 
in Visit 19. At the distance of the Coma Cluster, nearly 
all globular clusters are unresolved with HST, and can be 
well- approximated by point sources (the mean half-light 
radius of GCs, r^ « 3 pc, is only ~ 6% the full width 
half maximum of the point spread function). Because 
detection is done only in the F814W band, we only add 
artificial stars to these images. 

When adding point sources into the images, we avoid 
objects in the image as well as artificial stars already 
placed so as to avoid incompleteness due to confusion. 
We run the exact same detection pipeline on these im- 
ages as we do to create our object catalog and record 
whether the objects were detected as a function of mag- 
nitude and position. The number of artificial stars added 
and measured — 7,000,000 in Visit 19 alone, and approx- 
imately 4240 arcmin -2 for the other visits — ensures that 
we can derive a completeness curve for any position in 
the survey, and for any GC selection criteria. 

In a typical blank area in our images observed with 
the full exposure time, the 90% completeness level is at 
/ rj 26.8 mag, and the 50% completeness level is at J rj 
27.3 mag. At R ss 2' from the center of NGC 4874, 
however, these limits are 1.5 mag shallower. 

2.4. Globular Cluster Candidate Selection 

One of the main benefits of GC studies with HST is the 
ability to use morphology and resolution to separate GCs 
from their main contaminants, background galaxies. At 
Coma distances, GCs are point sources when observed 
with HST, but the great majority of background galax- 
ies are resolved. We use this ability to select against 
background contaminants and produce a relatively clean 
sample of GC candidates. 

We use a rough but effective concentration crite- 
rion to select GCs. Figure [2] shows the "magnitude- 
concentration" diagram for objects, where we measure a 
concentration index, 6*4-10 using the difference in mag- 
nitude measured in a 4 pixel diameter aperture and a 
10 pixel diameter aperture. Figure [2] shows that this 
index works well to distinguish point sources from ex- 
tended sources. Here, we show the distribution of ob- 
jects in Visit 19 (the one containing NGC 4874), which 
has the largest number of GC candidates. We overplot 
the objects from Visit 59, which is the most remote of 
our fields, and contains mostly background galaxies. The 
red lines show our selection region where we exclude 
nearly all of the background galaxies. Although we ex- 
perimented with different cuts in this diagram, includ- 
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Fig. 2. — / magnitude versus concentration index (C4_io = 
mipix — miOpix) for Visit 19, the central pointing containing 
NGC 4874 (left) and Visit 59, the background pointing most dis- 
tant from the cluster center which contains mostly background 
galaxies (right). The vertical locus of points around C4—10 = 0.45 
contains point sources. Most of the point sources in Visit 19 are 
likely to be GCs. The background galaxies are mostly resolved to 
be more extended than the GCs until I ~ 27, where some overlap 
the stellar locus. The red outlines shows our selection region for 
GCs. 

ing a variable width of the selection region with mag- 
nitude, the variations were not significant and we de- 
cided in the end that simplicity was best, choosing a cut 
of ±0.2 mag around the median concentration for point 
sources ((€4,-10) = 0.45). 

For the purposes of this study, we wish to maximize 
the number of good GC candidates, while also balanc- 
ing the increasing number of background contaminants 
with magnitude. Because of the depth and high spa- 
tial resolution of our data, we chose a fairly conservative 
magnitude limit, including objects with / < 26.5 mag. 
At this magnitude, our data is 97% complete in re- 
gions free of galaxy light, so completeness corrections 
are only important toward the centers of galaxies. At 
the Coma Cluster distance (m — M = 35), assuming 
an extinction Ai = 0.017 mag for NGC 4874 (Schlcgcl, 
Finkbeiner, & Davis 1998), this limit should include a 
significant fraction of the GCs in a Gaussian GC lumi- 
nosity function (GCLF) typical of giant ellipticals. We 
use the recently measured /-band GCLF measurement 
for the Virgo cD galaxy M87 (Peng et al. 2009), which 
was performed with deep HST/ ACS observations in the 
same F814W filter used by the ACS Coma Survey. Peng 
et al. (2009) quote a GCLF Gaussian mean and sigma 
of /-ii^vega = —8.56 mag and a = 1.37 mag. For AB 
magnitudes, we add 0.436 mag to fiiy e ga (Sirianni et al. 
2005). Assuming these values for a Gaussian GCLF, our 
GC catalog magnitude limit should include ~ 39% of all 
GCs, and ~ 75% of the luminosity in GCs. 

This is likely an oversimplification, however, as both 
the mean and width of the GCLF is known to vary with 
galaxy mass (Jordan et al. 2006; 2007). If we assume 
a Gaussian GCLF typical of dwarf ellipticals in clusters 
([ii.Vega = —8.1 mag and a = 1.1 mag, Miller & Lotz 
2007), then our limit includes ~ 22% of the total num- 
ber. This discrepancy is one of the main systematic un- 
certainties in our analysis. We emphasize, however, that 
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Fig. 3. — Spatial distribution of ACS GC candidates shown 
on a 1° X 1° Digitized Sky Survey image of the Coma Cluster 
with North up and East to the left (1.75 x 1.75 Mpc at Coma 
distance). At the top left of the image is the observed portions of 
the cluster core central mosaic. The largest concentration of GCs 
is around the central galaxy, NGC 4874. The other large galaxy, 
unobserved by ACS, is NGC 4889. At the bottom of the image 
are six fields in the outer regions of the cluster. The three outer 
fields to the right (west) show higher numbers of GC candidates 
because of their proximity to large galaxies. The eastern three 
outer fields (bottom center) are not near large galaxies and are used 
as background fields. The density of GC candidates throughout the 
entire cluster core is much higher than in the background regions, 
implying a large population of intracluster GCs. The ACS field 
sizes are roughly 202"on a side. 

changing the assumed GCLF does not affect the signif- 
icance of our result, just the inferred total number of 
GCs. Given that the depth of our data is not sufficient 
to measure the GCLF parameters directly, we choose to 
assume the brighter GCLF, seen in giant ellipticals, as 
this will give us a lower estimate for the number of GCs in 
any given area. The numbers could be higher by ~ 80% 
in regions where the GCLF for dwarf ellipticals is more 
representative. 

We also introduce a broad color cut of 0.6 < (g — I) < 
1.5 that should include all old globular clusters. This 
color range is based on the transformed g—z colors of 
GCs in the ACS Virgo Cluster Survey (Cote et al. 2004; 
Peng et al. 2006), and mainly eliminates distant, com- 
pact red galaxies. The ages of extragalactic GCs across 
all metallicities are primarily old (> 5 Gyr), especially 
those associated with massive early-type galaxies (e.g. 
Peng et al. 2004; Puzia et al. 2005; Beasley et al. 2008; 
Woodley et al. 2010), so this color range should include 
all bona fide GCs. 

3. SPATIAL DISTRIBUTION OF GC CANDIDATES 

Figure [3] plots the locations of GC candidates in our 
ACS images on a DSS image of the cluster. While it is not 
surprising that the number of GCs is high around massive 
ellipticals such as NGC 4874, what is striking about this 
figure is that the number of GCs across the entire central 
mosaic is high and is significantly elevated when com- 
pared to the numbers in the outer fields. Even the cor- 
ner fields of the central mosaic have many more GCs. Of 
the six outer fields, three in the southwest (lower right) 
have visibly elevated GC numbers due to their proximity 



Fig. 4. — Smoothed spatial distribution of GCs in the Coma 
Cluster core (30.'8 x 23.'0, 900 x 670 kpc). Pixels are 20" on a 
side, and color represents the surface density of GCs, corrected 
for completeness (blue to red denotes low to high density). The 
entire image has been smoothed by a Gaussian kernel with a = 
30". The dominant concentration of GCs is around NGC 4874, 
and an extended structure of GCs appears to connect NGC 4874, 
4889, and 4908. Some peaks in the distribution represent individual 
cluster galaxies. 

to NGC 4839 (top) and NGC 4827, two giant early-type 
galaxies. The three other fields to the south are not near 
massive cluster members. We take these three southern 
fields as an upper limit on the background contamination 
from foreground stars and compact distant galaxies. All 
of these fields have fewer GC candidates than does any 
field in the central mosaic. 

Other than an obvious concentration around 
NGC 4874 and NGC 4889 (the latter of which was 
not observed with ACS), the GC distribution is rela- 
tively uniform across most of the central mosaic and 
not spatially clustered; i.e., with the exception of the 
two central ellipticals, the spatial structure of the GCs 
is not highly correlated with the positions of cluster 
galaxies. This is partly a bias introduced by the failure 
of Source Extractor to detect GCs that are immediately 
in the vicinity of bright galaxies. However, GC detection 
should not be a problem in the halos of the galaxies, 
and except in a few cases we do not detect the kind of 
small-scale substructure one would expect in the cluster 
GC distribution if all the GCs were tightly associated 
with galaxies. 

Figure 2] shows more clearly the distribution of GCs 
in the cluster core. To produce this figure, we divide the 
core region into 20" x 20" "pixels" with each representing 
the surface density of GCs, corrected for spatially varying 
completeness, and smoothed with a Gaussian kernel with 
g = 30". The large concentration of GCs at the center- 
right is the GC system of NGC 4874. The GC system 
of NGC 4889 is also evident, although the galaxy itself 
was not observed. While Figure [3] shows that the overall 
surface density of GCs is well above the background, Fig- 
ure 2] shows hints of large-scale substructure in the GC 
spatial distribution. There appears to be an extended 
structure of IGCs connecting NGC 4874 to NGC 4889 to 
NGC 4908 and IC 4051, both of which lie just beyond 
the eastern edge of the mosaic. 

These observations suggest the existence of a large in- 
tergalactic population of globular clusters. In the follow- 
ing sections, we seek to verify and quantify their exis- 
tence. 
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4. BACKGROUND ESTIMATION AND GALAXY MASKING 

Contaminants to our sample of GCs consists of fore- 
ground stars and faint, unresolved background galax- 
ies (the sum of which we generically refer to as "back- 
ground" ) . This background is important to quantify, as 
a smooth background can mimic a smooth ICC popula- 
tion. Ground-based IGC studies in Coma are typically 
plagued by high background due to their inability to dis- 
tinguish distant galaxies from point sources. 

As a measure of our background, we choose the three 
outer ACS fields — visits 45, 46, and 59 — that are not 
near giant galaxies, and are shown at the bottom-center 
of Figure [3j For each of these fields, we select GC can- 
didates as described earlier, and also mask the regions 
containing a few obvious Coma members using the pre- 
scriptions described below. The surface density of GC 
candidates over these three fields is 2.8 ± 0.3 arcmin -2 , 
or ~ 28 per ACS field. As we will show, this is nearly an 
order of magnitude lower than the density of GCs even 
in the outer fields of the Coma core. 

To verify this background level, we compared our point 
source counts in these fields to those in the COSMOS 
HST Treasury project (Scoville et al. 2007). The COS- 
MOS survey imaged 1.8 deg 2 at high Galactic latitude 
with the same camera (ACS/WFC), filter (F814W), and 
depth as the Coma Cluster Survey. The number of point 
sources in our three background fields is entirely consis- 
tent with the numbers expected from the surface density 
of stars in the COSMOS fields. Down to Isu < 25 mag, 
we detect 104 ± 10 point sources in our three background 
fields and 98 are expected using the average surface den- 
sity from the COSMOS data. This independent check 
gives us more confidence that our background value is 
correct. 

The fact that the global background is so low com- 
pared to the detections in our Coma core fields gives us 
confidence that we are indeed detecting GCs within the 
Coma cluster. The more difficult question is whether 
these objects are truly "intergalactic" or simply part of 
extended galactic systems. This debate is not one easily 
resolved by imaging data alone. We can, however, ad- 
dress the contribution from galactic GC systems in two 
ways. First, we aggressively mask regions around known 
bright galaxies. Second, we can make certain assump- 
tions about the numbers and spatial extent of the GC 
systems of observed cluster members and compare sim- 
ulated GC distributions to the observations. We do this 
in order to test the hypothesis that the GCs observed in 
the cluster core are an intergalactic population. 

The details of these two methods are described in Ap- 
pendix [A] In short, we generate masks around all galax- 
ies with luminosities down to M g < —17 mag, both in 
and around all of our fields. The detection algorithm 
that we use for GCs actually ends up masking GCs 
around fainter galaxies because our chosen background 
estimation parameters cannot follow the steeply rising 
surface brightness profiles at the centers of galaxies. For 
each galaxy, we apply a liberal, size-dependent mask to 
the surrounding regions. These masks should eliminate 
~ 90% of the "galactic" GCs from our catalogs. For the 
remaining outer GCs, we subtracted a model GC system 
using an assumed Sersic n = 2 profile for the GC surface 
density, a reasonable assumption given previous measure- 



ments of GC system radial profiles. The parameters of 
this model are estimated based on scaling relations for 
S N and R e from Peng et al. (2008, 2011). This model- 
ing is done for all Coma galaxies in the Eisenhardt et al. 
(2007) catalog, which is complete to My < — 16 mag and 
extends to My < —14 m ag. This is described in greater 
detail in Appendix IA.2I 

Although we have taken great pains to model and sub- 
tract any residual GCs that may belong to Coma galax- 
ies, we find that our final result is largely insensitive to 
the assumed parameters. The detection of IGCs, as we 
will show below, is highly significant and not dependent 
on the details of the background or the modeling of GC 
systems. We estimate the systematic uncertainty due to 
our modeling procedure to be j^ooo GCs fAppendix lA.3p . 
only 5-9% of the inferred IGC population. 

5. RESULTS 

5.1. Radial Profile of GCs in the Coma Cluster 

NGC 4874 has previously been observed to have a large 
number of globular clusters (Blakeslee & Tonry 1995; 
Harris et al. 2000, 2009). It has also been shown to have 
a GC system whose spatial profile is shallower and more 
extended than those for other elliptical galaxies (Harris 
et al. 2009). Could the GCs that we see filling the cluster 
core simply be the extended GC system of NGC 4874? 
The situation is complicated by the fact that the core of 
the cluster contains not one but two giant ellipticals, the 
other being NGC 4889, as well as many other member 
galaxies. 

In Figure [5j we show the radial distribution of GCs 
in the cluster core, centered on NGC 4874. For each 
bin in radius, we sum up the number of observed GC 
candidates in unmasked regions, subtract the expected 
contribution of GCs from other Coma galaxies (shown 
as the dotted line in Figured]), and subtract the global 
background level (dot-dash line), leaving what should be 
the NGC 4874 and IGC population. We determine the 
mean completeness of the sample within the annulus, and 
extrapolate the total number of GCs assuming the M87 
GCLF as described above. We then sum the total ob- 
served, unmasked area within the annulus to determine 
the surface number density of GCs. The random errors 
in each bin are derived from the Poisson errors for the 
number of candidate GCs as well as from the Poisson 
error in the background, added in quadrature. 

This profile, tabulated in Table [TJ represents our best 
estimate of the radial surface density distribution of the 
GC system surrounding NGC 4874, uncontaminated by 
the GC systems of other cluster members. Perhaps the 
most interesting feature of this profile is a marked inflec- 
tion at R ~ 200 kpc, beyond which the GC surface num- 
ber density decreases much more slowly with radius. We 
interpret this flattening of the profile as the region where 
a large and extended population of IGCs starts to dom- 
inate the GCs directly associated with NGC 4874. The 
significance of this detection is extremely high, as the 
background level is shown in Figure [5] by the horizontal 
dot-dashed line at the bottom, with the estimated error 
of the background denoted as the shaded gray region. 
The inferred IGC surface density is a factor ~ 7 over the 
background. Another point of comparison is with the 
modeled surface density of remaining unmasked galactic 




Fig. 5. — The radial distribution of GCs in the Coma Cluster 
core centered on NGC 4874. The surface density of GCs in each 
bin (black points) is calculated after masking around known galax- 
ies, and statistical subtraction of GCs belonging to these cluster 
members. The radial profile exhibits a flat inner core as well as an 
inflection and flattening at large radii. We interpret the flat dis- 
tribution at large radii as evidence of a large population of IGCs. 
The dot-dashed line and gray band at bottom denote the surface 
density of background objects (plus la errors) determined from our 
outer ACS fields and subtracted from all radial bins. The back- 
ground level is a factor ~ 7 below that in the outermost bins. The 
arrow at bottom right shows the mean distance from NGC 4874 
of the three background fields. The dotted line is the modeled 
radial distribution of GCs belonging to cluster members that are 
still visible after masking. These have been subtracted from the 
GC radial profile, although they too are well below the overall level 
by a factor of a few. The data are well fit by a Sersic model plus 
a constant level (solid line). The Sersic component alone is shown 
as the dashed line. 

GCs, shown as the dotted line, which also has already 
been subtracted from our total GC profile. The overall 
surface density of GCs in this cluster profile is well above 
the surface density of masked galactic GCs (by a factor 
4-7), and thus the GCs we see are likely to be truly in- 
tergalactic. We have also found that our results do not 
change significantly if we only use data from the eastern 
or western half of the Coma core. 

A single Sersic profile, normally a good fit to the sur- 
face density profiles of GC systems, is not sufficient to 
describe the data for the central Coma Cluster GCs. In- 
stead, we fit a model combining a Sersic profile and a con- 
stant. It is likely that the IGCs have a radially decreasing 
density profile (although the simulations of Bekki & Ya- 
hagi (2006) suggest that they can also have a flat density 
distribution within the central few hundred kiloparsecs), 
but the data only allow us to measure their mean sur- 
face density. The solid line in Figure [5] traces the best 
fit model, and the dashed line that follows it until large 
radii is the best fit Sersic component. The fitted surface 
density of IGCs is 0.055 ± 0.002 kpc" 2 , which is a 19a 
detection over the background, 0.00845 ± 0.001 kpc -2 , 
assuming Poisson random errors. As we discuss in Sec- 
tion 2] and Appendix IA.31 we also need to account for 
systematic uncertainties from our modeling, but they do 
not affect the main conclusion, which is high significance 
of the detection. We list the best fit parameters in Ta- 
ble n 



TABLE 1 

GC RADIAL SURFACE DENSITY PROFILE 
CENTERED ON NGC 4874 



(R) 


<7 


Area 


(arcmin) 


(arcmin -2 ) 


(arcmin 2 ) 


0.036 


3276 ± 2321 


0.00546 


0.058 


4031 ± 2331 


0.00719 


0.082 


1856 ± 1076 


0.01428 


0.115 


2175 ± 825 


0.02834 


0.162 


2950 ± 507 


0.05617 


0.229 


2737 ± 326 


0.11152 


0.322 


1978 ± 181 


0.22114 


0.454 


i 0Q7 1 inn 

183/ ± 128 


0.38152 


0.639 


1423 ± 91 


0.56911 


0.900 


1164 ± 71 


0.73770 


1.268 


828 ± 52 


0.96299 


1.786 


641 ± 39 


1.45523 


2.516 


318 ± 15 


4.62493 


3.543 


161 ± 7 


12.59246 


4.991 


97 ±4 


22.04740 


7.030 


62 ± 4 


16.82656 


10.066 


44 ±3 


19.62118 


13.674 


53 ±2 


40.98172 


17.886 


38 ±3 


14.17776 


Note. 


— These surface densi- 



ties are corrected for completeness, the 
full Gaussian GCLF, and include the 
masking and subtracting of GCs be- 
longing to other cluster members, as 
described in Section 15.11 The surface 
density of contaminants (also correct- 
ing for the GCLF) as marked by the 
dot-dashed line in Figure[5]is 7.2 ±1.2. 

Although we cannot determine the shape of the IGC 
component's density profile, one constraint is that it 
must fall rapidly after the limits of our data. The mean 
distance of the three fields we are using to measure the 
background is shown as the vertical arrow at 1.3 Mpc. 
Therefore, the GC surface density must fall to zero, or 
at least the level of the dashed line, by this distance. A 
steep falloff like this favors a low n Sersic profile for the 
IGCs (n = 1-2), similar to the ICL profiles in Seigar et al. 
(2007) and X-ray gas in galaxy clusters (Demarco et al. 
2003), but lower (i.e., steeper in the outer regions) than 
dark matter halo density profiles (Merritt et al. 2006). 
However, at these radii, it may not make as much sense 
to speak of a circularly symmetric GC radial profile, and 
it would be more useful to map in two dimensions the 
spatial distribution of GCs. 

5.2. Total Numbers of GCs and Specific Frequency 

We use our radial spatial density profile to estimate the 
total number of GCs in the cluster core, which wc define 
to be the extent of our data. Integrating this profile 
for R < 520 kpc gives a remarkable 70000 ± 1300 GCs, 
with the IGC component dominating the GC population 
beyond 150 kpc. As listed in Table [2j the number of 
GCs belonging to NGC 4874's "Sersic component" out to 
this radius is « 23000, leaving a remaining 47000 ± 1600 
(random) ^5000 (systematic) to be IGCs. There are over 
twice as many IGCs as there are GCs from the Sersic 
component, resulting in an IGC fraction of the entire 
central GC system of ~ 70%. 

With a measurement of NGC 4874's luminosity, we can 
calculate an "intrinsic" specific frequency for the galaxy. 
Harris et al. (2009) use a luminosity of My — —23.46 
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TABLE 2 

Best parameters for Sersic plus constant (£k?c) model 
fit to Coma central GC system 



Parameter 


Value 


Description 


n 


1.3 ±0.1 


Sersic index 


Re 


62 ± 2 kpc 


Sersic effective radius 


Se 


0.437 ±0.034 kpc" 2 


GC surface density at R e 




0.055 ± 0.002 kpc" 2 


Mean IGC surface density 


NGC,tot 


70000 ± 1300 


Total GCs within 520 kpc 


NgC, Sersic 


23000 ± 700 


"Sersic" GCs within 520 kpc 


NlGC 


47000 ± 1600 (r) +_fZ (s) 


IGCs within 520 kpc 



(adjusted to D — 100 Mpc), but surface brightness pro- 
files from SDSS imaging (J. Lucey, private communica- 
tion) and KPNO 4-meter CCD imaging (R. Marzke, pri- 
vate communication) show the galaxy to be substantially 
brighter. Measurements of the total r-band luminosity 
from mosaicked SDSS frames gives a value of r = 10.23 
(M r = -24.77, assuming E(B-V) = 0.009 from Schlegel 
et al. 1998), which includes a 0.32 mag extrapolation us- 
ing the best-fit Sersic profile for the light beyond R = T . 
The mean color of the galaxy is g — r ss 0.8 mag, which 
produces My — —24A7 using the Lupton (2005) trans- 
formation. 

With this, we calculate an "intrinsic" specific fre- 
quency for NGC 4874 of S N = 3.7 ± 0.1 (the errors are 
purely from the total numbers of GCs and do not include 
errors in the luminosity). 2 This value is very much in 
line with the those of non-cD, giant early- type galaxies 
in the Virgo and Fornax clusters. 

If we assume that all GCs (including IGCs) are part 
of the NGC 4874 system and that we are not missing 
any luminosity from the galaxy, then this would give the 
specific frequency within 520 kpc a higher value of Sn = 
11.4 ± 0.2, a value similar to those measured for some cD 
galaxies. Another interpretation, which we discuss later, 
is that the specific frequency of the system is the lower, 
more normal value, but that the IGCs are tracing a large 
amount of intracluster star light that is unaccounted for. 

5.3. Comparison to NGC 4874 Surface Brightness 
Profile 

A relevant comparison for the GCs is to the surface 
brightness profile of the field star light of NGC 4874. In 
Figure [6j we plot the light profile in circular apertures 
around NGC 4874 from two independent data sets (with 
arbitrary normalization). As mentioned in the previous 
section, the first is from measurements using SDSS r- 
band imaging (J. Lucey, in preparation), and the second 
uses imaging from the Mosaic-I camera on the KPNO 
4-meter telescope (R. Marzke, private communication). 
Both profiles are in good agreement in the inner regions 
(R < 20 kpc), but start to diverge in the outer regions 
due to differences in the sky measurements. The differ- 
ence between the two profiles is at the level of 2% of the 
sky. 

In the regions beyond 100 kpc, we also plot the surface 

2 If we use the older, fainter value for the luminosity, then 
Sjv = 9.5 ± 0.3, which is consistent with the value found in the 
HST/WFPC2 study of Harris et al. (2009). Their data only ex- 
tended to R ~ 65 kpc, and thus were not able to detect the IGC 
population. The higher value is also consistent with the value esti- 
mated by Blakeslee et al. (1997), who also used a fainter luminosity. 
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Fig. 6. — The radial distribution of GCs centered on NGC 4874 
(black dots) compared to the surface brightness profile of field star 
light around NGC 4874. Three different sources arc used for the 
surface brightness profiles: KPNO (solid), SDSS (dotted), Thuan 
& Kormendy (1977, dashed). The shaded and striped regions rep- 
resent a change in sky determination of ±2% for the KPNO and 
SDSS profiles, respectively. The GC radial surface number density 
profile and the field star surface brightness profile do not exhibit 
similar shapes at either small or large radii. The surface brightness 
measurements at large radii, however, are entirely dependent on an 
accurate measure of the sky brightness. 



brightness profile of the "intracluster background light" 
in the Coma cluster as determined photographically by 
Thuan & Kormendy (1977), transforming from Q to r 
magnitudes using an offset of (Q — r) = 0.37 mag, based 
on a (B — V) = 0.7 mag, their published transformation 
from Q to Thuan-Gunn r (Thuan & Gunn 1976), and 
then an offset to SDSS r (Fukugita, Shimasaku, Ichikawa 
1995). This profile appears to match the SDSS photom- 
etry at R ~ 100 kpc but then continues with a shallower 
slope. 

It is clear that at large radii (R > 100 kpc), the de- 
termination of the sky is crucial to the measurement of 
the ICL. To illustrate this, we shade in the regions cor- 
responding to a change in sky determination of ±2% of 
the sky level around the KPNO and SDSS measured pro- 
files. Although there is no evidence for a "break" in the 
measured surface brightness profiles akin to what we see 
in the GCs, the surface brightness profile in the regions 
where IGCs dominate GC counts is entirely dependent 
on the determination of the sky level to better than 1% 
and thus is difficult to quantify. 

We can use these profiles to calculate the "local" spe- 
cific frequency of the outer GCs, although any calcula- 
tion is highly uncertain due to sky subtraction for the 
surface photometry. Nevertheless, we can take these sur- 
face brightness profiles at face value to see if the cal- 
culated values are reasonable. Assuming that the sur- 
face brightness at R = 200 kpc is /i r ~ 26.2 mag 
(following the Thuan & Kormendy profile), and using 
V — r = 0.2 mag for old metal-poor stellar populations, 
then /iy(200 kpc) ~ 26.4 mag. Given the IGC sur- 
face density at these radii (46 arcmin -2 ), we estimate 
5^(200 kpc) = 5. If we assume the profile derived from 
SDSS data, however, then ^y(200 kpc) w 27.5 mag, re- 
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suiting in Sn (200 kpc) = 13. We emphasize that the 
surface photometry is extremely uncertain at these radii, 
and the upper error bar on this number is essentially 
unconstrained. The local values of Sn at these kinds 
of radii have previously been reported to be quite high 
(Tamura et al. 2006 in M87 and Rhode & Zepf 2001 for 
NGC 4486), but those measurements are equally uncer- 
tain for similar reasons. 

In the inner regions, there is a notable divergence be- 
tween the GCs and the galaxy light. The galaxy does 
not show the prominent core within 10 kpc that the GC 
system does, only displaying a flattening in the profile at 
a smaller radius. 

5.4. Comparison to the M87 GC System 

Perhaps the most relevant local comparison for the 
NGC 4874/Coma Cluster GC system is that of M87 in 
the Virgo Cluster. In Figure [7] we show the GC radial 
surface density profiles of the two GC systems. The 
outer M87 profile is taken from the data of McLaugh- 
lin (1999) and Tamura et al. (2006), while the central 
regions of the profile are from the ACS Virgo Cluster- 
Survey (ACSVCS) data shown in Peng et al. (2008). We 
note that the physical resolution of the Coma .ff ST data is 
very competitive with ground-based Virgo observations 
(O'.'l resolution at Coma distance is equivalent to 0'.'6 res- 
olution at Virgo distance), but obviously cannot match 
the ACSVCS observations of M87. None of the Virgo 
data sets go as far out in physical radius as our Coma 
data, but they still provide a useful comparison. The two 
GC systems profiles are similar in the range of intermedi- 
ate radii (20-100 kpc), but differences appear in the very 
inner and outer regions. Most noticeably, the Coma GC 
systems displays a very pronounced core within 10 kpc, 
which does not appear to be present in the Virgo sys- 
tem except perhaps within 1 kpc. This deficit of GCs 
at the center of NGC 4874 is also evident in the analy- 
sis of Harris et al. (2009). This core could be the result 
of dynamical friction destroying GCs at the center of 
NGC 4874. 

The core in the GC profile, and the divergence from the 
M87 GC profile, is most evident within 10 kpc. Could 
this be due to unaccounted observational incomplete- 
ness? The four innermost radial bins (R < 3.5 kpc), have 
the largest errors and shallowest observations because of 
the bright galaxy light (completeness of the GCLF is 
~ 25% in these bins). There are multiple reasons, how- 
ever, why we believe these lower surface densities to be 
real. The radius at which the core becomes apparent, 
20", is large for HST imaging. Even excluding the inner 
10" (4.8 kpc), the difference in slope between the two GC 
profiles is still apparent. Also, the completeness tests we 
apply also take into account incompleteness due to im- 
perfect profile subtraction, and thus is a true measure of 
the completeness in these radial annuli. In order to turn 
the M87 GC profile into the Coma GC profile through a 
systematic overestimation of the completeness, the com- 
pleteness would have to be overestimated by nearly an 
order of magnitude. Lastly, this deficit of GCs relative to 
the galaxy light profile was also independently found by 
Harris et al. (2009) using HST/WFPC2 data (see their 
Figure 6). 

In the outer regions the M87 profile follows the single- 
Sersic fit out to the limits of the data. The Virgo data, 
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Fig. 7. — The radial distribution of GCs centered on NGC 4874 
(black dots) compared to the distribution of GCs around the Virgo 
cD galaxy M87. M87 data is from McLaughlin (1999, McL99, 
orange diamonds). Tamura et al. (2006, T06, red asterisks), and 
the ACS Virgo Cluster Survey (Peng et al. 2008; blue triangles) 
with a Sersic fit to the combined data set overplotted (dot-dashed). 
The Coma GCs have a much shallower and larger core, as well 
as an inflection where IGCs begin to dominate. Even with the 
larger T06 data set, there is not yet evidence for a profile inflection 
around M87 like what we see in Coma, although the data do not 
go comparably far out in radii. 

however, only reaches a radius of 130 kpc, and thus would 
not be sensitive to the kind of IGC population we see in 
Coma. In fact, if the Coma data had the same physical 
radial extent, it would have been very difficult to de- 
tect the IGC population. More deep, wide-field imaging 
of the area around M87, such as the Next Generation 
Virgo Survey 3 , will be necessary to detect or place more 
stringent limits on a population of IGCs in Virgo. 

5.5. GC Color Distributions 

For old GCs, the broadband color is an indicator of 
metallicity. The color distributions of extragalactic GC 
systems have been studied extensively with HST (e.g., 
Larsen et al. 2001; Peng et al. 2006), and are often bi- 
modal in nature (Gebhardt & Kissler-Patig 1999), espe- 
cially in massive early-type galaxies. In FigureJSJ we plot 
the color distribution of bright GCs (applying a magni- 
tude limit of g < 25 mag for higher S/N) in the unmasked 
regions of the Coma Cluster core. The color distribution 
of all GCs shows the typical bimodality seen in extra- 
galactic GC systems, displaying a prominent peak of blue 
(metal-poor) GCs with (g — /) sa 0.9, and a red (metal- 
rich) peak with (g — I) w 1.15. 

We plot the color distributions divided by distance 
from the center of NGC 4874 — those within 50 kpc 
(galactic GCs) and those outside of 130 kpc (predomi- 
nantly IGCs). We use the Kaye's Mixture Model (KMM; 
McLachlan & Basford 1988; Ashman, Bird, & Zepf 1994) 
implementation of the expectation-maximization (EM) 
method to fit two Gaussians with the same standard de- 
viation to the GC color distributions of each sample. 
Both the galactic and intergalactic GCs are much bet- 

3 The Next Generation Virgo Survey (NGVS) is a Large Program 
with the Canada-France-Hawaii Telescope. 
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Fig. 8. — The color distributions of all GC candidates in un- 
masked regions with I < 25 (black solid line), a predominantly 
IGC subsample with R > 130 kpc (blue dashed line), and a pre- 
dominantly galactic GC subsample with R < 50 kpc (red dotted 
line). The total distribution exhibits a bimodality typical for ex- 
tragalactic GC systems, as does the IGC sample. For the IGCs, 
blue GCs outnumber red GCs by a ratio of 4:1. Only in the inner 
regions, where the GC system of NGC 4874 is dominant, does the 
number of red GCs compare to the number of blue GCs. The color 
distribution of the inner GCs is also bimodal, but with the blue 
population having a much redder color. 



ter described by bimodal distributions in color than by a 
single Gaussian with p-values less than 0.001. The inner 
blue GCs, however, have a much redder mean color — 
(g - I) = 0.94 as opposed to (g - I) = 0.89 for GCs 
in the outer regions — such that they are nearly merged 
with the red GCs (g - I) = 1.18. the metal-poor GCs 
are either quite red or there is a substantial population 
of GCs at intermediate color. It could also be the sign 
of a significant radial color gradient within the blue GC 
subpopulation (e.g., Harris 2009, Liu et al. 2010). This 
gradient could be in metallicity or in age, reminiscent of 
younger metal-poor GCs in Local Volume dlrr galaxies 
(Sharina et al. 2005, Georgiev et al. 2008, 2009), although 
the colors of the blue GCs are not so blue as to require 
very young ages. The inner regions also have equal num- 
bers of blue and red GCs where the fraction of red GCs 
is fred = 0.51. 

The generally red colors of the inner GCs was noted 
by Harris et al. (2009), who used their WFPC2 data to 
show that the inner regions of the GC system had a high 
(> 50%) fraction of red GCs. Using the larger ACS cov- 
erage of our survey, however, we see that the total GC 
system around NGC 4874 is dominated by blue GCs. In 
the outer sample plotted in Figure the IGC popula- 
tion is dominated by blue GCs with a blue-to-red ratio 
of 4:1, (fred — 0.2) although the fraction of red GCs is 
still significant. Even measuring only the Sersic compo- 
nent of the GC subpopulations, the NGC 4874 GCs are 
dominated by blue GCs, with the red GCs being more 
spatially concentrated. The radius at which the surface 
density of blue and red GCs are equal is R « 40 kpc. 
Thus, like the giant ellipticals in the Virgo Cluster (Peng 
et al. 2008), the GC system of NGC 4874 as well as the 
IGC population, is dominated by blue GCs. 



The mean colors of the blue and red outer GCs are 
(g — I) = 0.89 and 1.22, respectively. The blue peak 
has a typical color for the GC systems of galaxies with 
masses at or below L* (Peng et al. 2006). The red peak, 
however, is still fairly red, having nearly the same mean 
color as the inner GCs around NGC 4874. We caution 
that some of these red GCs may be from the large metal- 
rich GC system of NGC 4889 and other more massive 
cluster members. 

We can do a simple test to see if the red GCs beyond 
130 kpc are associated with luminous galaxies. We calcu- 
late the distance from each GC to the nearest luminous 
(Mb < —17) galaxy, and compare the mean shortest dis- 
tance for the red and blue GCs. We find that there is 
no significant difference between these two populations. 
The median shortest distance for red GCs is 75" and that 
for blue GCs is 74", with the biweight Gaussian sigmas 
of each distribution being 35". The spatial behavior of 
the red and blue GCs beynd 130 kpc are identical and is 
consistent with a population of IGCs mostly uncorrelated 
with nearby galaxies. 

6. DISCUSSION 

6.1. Intracluster Light Inferred From IGCs 

The existence of intracluster starlight and globular 
clusters has for decades been considered important, but 
always difficult to observe. Recently, as observations 
have improved and theory has shown that intracluster 
stellar populations are an essential feature of galaxy- 
galaxy and galaxy-cluster tidal interactions, there is in- 
creased interest in quantifying its properties: total mass, 
spatial distribution, metallicity, and kinematics. We 
have shown that with HST, a direct detection of IGCs 
is a clean way to measure one component of intracluster 
stellar populations in nearby galaxy clusters. 

Only a small fraction of the total or stellar mass is in 
the form of old globular clusters. The total mass fraction 
in GCs appears to be relatively constant across galaxies 
and galaxy clusters (McLaughlin 1999; Blakeslee 1999). 
Because of the variation in the stellar mass-to-light ratio 
across galaxy mass, however, the specific frequency (or 
stellar mass fraction) can vary widely. Massive ellipticals 
and dwarf galaxies can have the highest Sn, and galaxies 
with luminosities around L* have the lowest Sn- The 
stellar mass fraction in early-type galaxies ranges from 
~ 0.2% to a few percent (Peng et al. 2008). 

We do not know the fraction of stellar mass that is in 
the form of IGCs, but we can make reasonable assump- 
tions. One possibility is that the IGCs originate from 
low-mass dwarf galaxies that are tidally disrupted. The 
Sjv of such a population can vary, and can depend on the 
clustercentric radius (Peng et al. 2008), but if we assume 
that dwarfs in the cluster center will have high specific 
frequencies, then we can reasonably assume Snjgc = 8, 
which is similar to the Sn values for high-SV dEs at the 
center of the Virgo Cluster (Miller & Lotz 2007; Seth 
ct al. 2004; Peng et al. 2008), the GC systems of even 
lower mass dwarfs (e.g., Puzia & Sharina 2008), and for 
Virgo cluster IGCs (Williams et al. 2007). 

Such a value for the specific frequency of the IGC 
population would imply a total ICL luminosity within 
520 kpc of My = —24.4 mag, an amount of star light 
equal to the whole of NGC 4874. Spread out uniformly 
over this entire area, this luminosity would have a mean 
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TABLE 3 

ICL Luminosity, Surface Brightness, Mass, and Fraction 
(R < 520 kpc) as a Function of Sn,igg 
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surface brightness of fiv,iCL = 27 mag arcsec -2 , which is 
still challenging for surface photometry, but detectable in 
the best observations (Mihos et al. 2005). If we assume 
M/Ly — 1.7, a value typical of cluster dEs, then this 
corresponds to a stellar mass of Micl ~ 9 x 10 11 Mq. 

Assuming a lower specific frequency like Snjgc = 1-5, 
as is more common for L* early-type galaxies or low-Sjv 
dEs and dSOs in the outskirts of the Virgo Cluster, would 
imply a higher luminosity for the ICL. Conversely, as- 
suming a very high specific frequency like that in M87 
or the highest-S^r dwarfs, Snjgc — 12, would imply 
a lower ICL luminosity. The inferred ICL luminosities, 
surface brightnesses, masses, and central ICL fractions 
for these three assumed values of Snjgc ar e listed in 
Table [3] These values bracket the range of "local" Sn 
we calculate from the outer surface brightness profiles of 
NGC 4874 in Section 15.31 Deeper, systematically con- 
trolled photometry of the intracluster region will be the 
best way to set a meaningful limit on the Sn of the IGC 
population. 

We have also calculated these values using the total 
luminosity in GCs, which is more robust than the total 
number because most of the luminosity in GCs is brighter 
than the mean of the GCLF (Harris 1991). Using typical 
values for the stellar luminosity fraction (0.2-1.2%) from 
Peng et al. (2008), we arrive at very similar numbers to 
those in Table [3] 

All of these inferred luminosities and masses are larger 
than those previously inferred from direct measurement 
of the low surface brightness ICL in Coma. Both Gregg 
& West (1998) and Adami et al. (2005) estimate a value 
of Mr ~ —22, which for an old stellar population is 
roughly My ~ —21.5. However, their values were the 
result of combining the luminosities of individual ICL 
sources. The 1% detection limit of the Adami et al. 
(2005) study is quoted as fiy = 25.8, which is much shal- 
lower than what is required if Sn.igc ^ 2. It is likely 
that these previous studies were sensitive to overdense 
regions in the ICL, but not to the overall population of 
intracluster stars. 

6.2. The Intracluster Fraction 

When considering only the IGC and NGC 4874 GC 
system (equivalent to BCG+ICL measurements), we find 
that the IGCs make up ~ 70% of the total "central" GC 
system. Although more difficult, we can also estimate 
the IGC fraction for all GCs in the cluster core, includ- 
ing the GCs belonging to cluster member galaxies. We 
do this in two ways. First, we can calculate this number 
only for those areas that we have observed. We assume 
the IGC surface density to be the fitted value in Tabled 
We can subtract the surface density of IGCs from our 
data and what remains is the "galactic" GC population 



in the observed areas. Using this metric, we find that 
roughly ~ 30% of the observed GCs belong to the intr- 
acluster population. This number does not account for 
the GCs close to galaxies that were missed by our de- 
tection algorithm, so it can be considered an upper limit 
over the area observed. 

The main drawback of limiting to the observed area, 
however, is that we do not observe large portions of 
the cluster core, including the other massive elliptical, 
NGC 4889. Another method is to apply the GC sys- 
tem modeling previously described to estimate the total 
number of GCs within 520 kpc that belong to galaxies. 
When we do this, we find that the IGCs make up ~ 45% 
of the GCs in the Coma Cluster core (this fraction is 
higher because our observed area includes NGC 4874, 
depressing the IGC fraction). This number is naturally 
more uncertain because a large fraction of the galaxies 
are unobserved. 

Given the assumptions in the previous section, we can 
infer that approximately half of the stellar luminosity 
associated with the "central" system (NGC 4874 and the 
ICL) is in the ICL. Because the ICL is likely to be more 
metal-poor and have a lower mass-to-light ratio than the 
stars at the center of NGC 4874, this translates to about 
one-third of the stellar mass being in the ICL. We would 
expect that the IGC fraction should be higher than the 
ICL fraction, because most of the stellar mass is in L* 
galaxies that have low Sn, whereas the specific frequency 
of the IGC population is likely to be high. This is why 
we find that the IGC fraction of the "central" system is 
~ 70%, but that the likely ICL fractions in Table [3] are 
closer to 50%. 

All of these numbers imply a high ICL fraction in the 
cluster, but are consistent with expectations from obser- 
vations and simulations for the core of a massive, rich 
cluster such as Coma. Seigar et al. (2007) find an central 
ICL fraction of 60-80% around cD galaxies, and Purcell 
et al. (2007) expect similar ICL fractions from simula- 
tions of tidal stripping. 

Previous studies of extragalactic GC systems have also 
suggested that there is a roughly constant ratio between 
the total number or mass of GCs in a system and the 
total dynamical or baryonic mass (e.g., Blakeslee 1999; 
McLaughlin 1999; Peng et al. 2008; Spitler et al. 2008). 
If this is the case, than the IGC fraction may be more 
representative of the total mass fraction in the diffuse 
component of the galaxy cluster. 

6.3. Scenarios for the Origin of IGCs 

Both theory and simulations suggest that an intraclus- 
ter stellar component is a natural result of many physical 
process that shape galaxies in the cluster environment. 
Galaxy-galaxy tidal interactions (e.g., Gallagher & Os- 
triker 1972; Moore et al. 1996; Stanghellini, Gonzalez- 
Garcia, & Manchado 2006), tidal forces as galaxies orbit 
through the cluster gravitational potential (e.g., Mer- 
ritt 1984; Gnedin 2003), tidal "preprocessing" within 
infalling galaxy groups (Rudick et al. 2006), tidal de- 
struction of low mass galaxies (Lopez-Cruz et al. 1997), 
and tidal tails that escape from merging galaxies (e.g., 
Murante et al. 2007) are all plausible mechanisms to lib- 
erate stars and GCs from the gravitational potential of 
their host galaxy. The problem is in trying to distin- 
guish between these mechanisms, and determine what 
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combination produces the observed intracluster stellar 
populations. 

Based on previous ICL studies in Coma, the specific 
frequency of the ICC population is not likely to be very 
low otherwise its associated ICL would already have been 
detected at relatively high S /N. One straightforward pos- 
sibility is that the GCs originate mainly from disrupted 
dwarf galaxies, which can have high specific frequencies. 
Dwarf galaxies, having low masses and surface mass den- 
sities, are prone to be destroyed by interactions with 
other galaxies and with the cluster potential. If we as- 
sume Snjgc = 8, then the inferred ICC population 
would require the disruption of <~ 2300 dwarf galaxies 
with My = -16. 

Simulations suggest, however, that the ICL originates 
from stars in L* galaxies rather than from dwarfs (Bckki 
& Yahagi 2006; Purcell et al. 2007), and that a signif- 
icant fraction of the ICL is first liberated in dynami- 
cally cold streams possibly thrown off from merger rem- 
nants (Rudick et al. 2009). ICL production mechanisms 
that depend on galaxy-galaxy interactions appear to be 
able to produce the right amount of intracluster mass 
(early calculations by Gallagher & Ostriker (1972) of 
stars stripped from L* galaxies during close encounters 
predicted that the total ICL mass in the Coma cluster 
should be approximately 10 12 Mq, a number quite close 
to what we infer from the IGCs). 

Galaxy-galaxy harassment, however, is expected to 
produce a more centrally concentrated ICC /ICL compo- 
nent than, for instance, interactions from with the mean 
cluster tidal field (Merritt 1984), because the cross sec- 
tion for galaxy-galaxy interaction is strongly peaked at 
the cluster center. We unfortunately cannot determine 
the full spatial density profile of the IGCs as our obser- 
vations only reach R = 520 kpc. What we do measure in 
the IGC dominated regime (150 < R < 520 kpc), how- 
ever, is consistent with a flat profile and shows no sign 
of being very centrally concentrated. 

Early- type galaxies near L* have a nearly universally 
low specific frequency, with Sn ~ 2 (Peng et al. 2008). If 
these kinds of galaxies were entirely disrupted to supply 
the necessary IGCs, it would imply a relatively high ICL 
surface brightness. What is more likely is that it is pre- 
dominantly the halos of these galaxies that are stripped, 
and that the stellar halos have higher local Sn- A pre- 
diction of this scenario would be the presence of a large 
number of intermediate-luminosity galaxies in the clus- 
ter core that have relatively few GCs and compact GC 
systems. These galaxies might be detectable as outliers 
in the color-magnitude relation for early-type galaxies. 

The color distribution of the IGCs also provides a clue. 
Although the IGCs are mostly blue, 20% of the IGCs are 
still red. This appears contrary to the simulations of 
Bckki & Yahagi (2006), which predicted a single metal- 
poor peak for the GCs with few metal-rich IGCs. The 
red IGCs we observe have a mean color consistent with 
the red GC populations found in the more massive early- 
type galaxies rather than in the dwarfs (Larsen et al. 
2001; Peng et al. 2006). The fraction of red IGCs is also 
slightly higher than that in dwarfs. For Virgo early-type 
dwarfs with My > —18, the mean fraction of red GCs is 
only 10%, and many have zero red GCs. Producing this 
many red IGCs solely from the disruption of low mass 
galaxies would be difficult. Stripping of GCs from the 



halos of more massive galaxies (My ~ —19) might be a 
more natural mechanism for producing the IGC popu- 
lation, although the stripping could not be too efficient 
otherwise the ICL might be too red (c.f., the blue ICL 
color seen in Virgo by Rudick et al. 2010). 

Mergers of these galaxies could also cause some metal- 
rich GCs to escape into intracluster space , although we 
are not able to differentiate between merging and strip- 
ping. Dissipational merging, with star formation, could 
also produce IGCs within gaseous tidal tails (e.g., Knapp 
et al. 2006), although for this mechanism to be impor- 
tant it would have to happen early in the formation of 
the cluster. 

Although it is likely that dwarf galaxy disruption does 
play a role in supplying IGCs, the non-negligible frac- 
tion of red IGCs indicates that stripping and merging of 
intermediate-mass galaxies contributes a significant frac- 
tion of the population. Unfortunately, we do not cur- 
rently have sufficient data to derive the full spatial profile 
of the IGC population itself, which would help determine 
the originating galaxy population, differentiate between 
merging, galaxy-galaxy interactions, and galaxy-cluster 
interactions, and also test whether the IGC radial den- 
sity follows the overall mass density. More ACS imaging 
within the core and at larger radii would be extremely 
useful for this purpose. 

7. CONCLUSIONS 

We present the spatial distribution and color distribu- 
tion of the globular cluster population in the Coma Clus- 
ter core, from imaging obtained as part of the HST/ACS 
Coma Cluster Survey. 

• We discover a large population of GCs that do not 
appear to be associated with individual galaxies. 
After masking and statistically subtracting the GC 
systems of all member galaxies except the central 
galaxy, NGC 4874, we find that IGCs dominate the 
GC surface number density at galactocentric radii 
beyond R > 130 kpc. These IGCs appear to have a 
flat surface density profile out to the extent of our 
data (R = 520 kpc). 

• Using a Sersic plus constant model fit to the radial 
profile, we estimate that there are 47000 ± 1600 
(random) (systematic) IGCs within 520 kpc, 
and that these make up 70% of the central Coma 
Cluster GC system (NGC 4874+IGCs). Including 
the GC systems of cluster members, we still esti- 
mate that IGCs make up 30-45% of all GCs in the 
cluster core. 

• The color distribution of the IGCs is bimodal, with 
red GC fraction of 20%. The inner GCs around 
NGC 4874 also have a bimodal color distribution, 
although with a fairly red metal-poor population, 
and with equal numbers of blue and red GCs. 

• The reasonably high red fraction (20%) for the 
IGCs, and red colors of the metal-rich IGCs, sug- 
gests that the IGC population did not solely origi- 
nate in dwarf galaxies, but at least some part was 
stripped from the halos of more massive galaxies. 

• These IGCs trace a large population of stars 
with an estimated surface brightness of [iy « 
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27 mag arcsec -2 (assuming Snjgc = 8). The ICL 
makes up approximately half of the stellar luminos- 
ity and one-third of the stellar mass of the central 
system. The IGCs and ICL are associated with the 
build up of the central cluster galaxy, and are likely 
the result of the continued growth and evolution of 
the cluster well after star formation has ceased. 
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APPENDIX 

A. MASKING AND MODELING OF GLOBULAR CLUSTERS FROM COMA GALAXIES 

A.l. Generating Masks 

We first generate masks around bright galaxies in and around all of our fields. We define masked galaxies to be 
those with 200" of an ACS field center, and with go < 18 (M g < —17) as observed in the SDSS Data Release 6 
(DR6) catalog. To be conservative, we choose all galaxies that meet these criteria, although nearly all at these bright 
magnitudes are cluster members. We are required to use a catalog external to our ACS survey because some large 
galaxies not observed by ACS, but just off the edge of a field, can contribute GCs to our observations. 

For each galaxy, excluding the two giants NGC 4874 and 4889, we use as a crude measure of galaxy size the radius 
that encloses 50% of the Petrosian light (Rp.so). After experimenting with various mask sizes, we liberally mask area 
corresponding to a circular aperture with a radius of 8i? p ,5o- For a Sersic profile with n = 2, this would mask 99% of 
the galaxy light. Even if the GC systems were twice the size of the galaxies, this would still mask ~ 92% of the GC 
system. For three of the larger galaxies, the SDSS sizes are not reliable, so we manually adjusted the mask size. For 
NGC 4889, IC 4045, and 4908, we used mask radii of 230", 68", and 160", respectively. 

A. 2. Simulated GC spatial distributions from galaxies 

For NGC 4889 in particular, the GC system is so large that a simple masking will not eliminate the contribution 
of its GCs. This is also the case for a few other bright galaxies. There is also the possibility that smaller but more 
numerous GC systems from the many dwarf galaxies in the cluster are contributing to the observed distribution of 
GCs. To account for this, we have simulated the composite GC distribution we would expect from the known galaxies 
in Coma. 

We start with the Coma Cluster galaxy catalog of Eisenhardt et al. (2007), who published UBVRIzJHK s photom- 
etry for Coma galaxies covering a region well-matched to our cluster core mosaic. For this purpose, their catalog is 
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Fig. Al. — Smoothed spatial distribution of GCs in the Coma Cluster core as in Figure [4] (30.'8 X 23f0, 900 X 670 kpc), except with 
areas around bright galaxies masked, and modeled contribution from galactic GC systems subtracted. The area and surface densities in 
this image are used in the radial profile shown in Figure [5] 

superior to the SDSS data because Eisenhardt et al. (2007) specifically measured photometry of the Coma galaxies, 
many of which are too large for the SDSS pipeline to handle properly. Their catalog is complete to My < — 16 mag and 
extends to My < —14 mag. We use their Tables 1 and 2 as the basis for our simulated GC systems. We assume that 
each galaxy in their catalog has a GC system that is circularly symmetric, and has a radial number density distribution 
that follows a Sersic (1968) profile with index n = 2, which is a good fit to the GC spatial density profiles of the Virgo 
giant ellipticals M87 and M49. The one exception to this is NGC 4889, for which it seems clear from Figure [3] that 
its associated GC systems is not spherical. For this galaxy, we assume that the spatial density distribution follows a 
Sersic profile along the geometric mean radius, (r) = \fab~ where a and b are the major and minor axes, respectively. 
The ellipticity of NGC 4889 is e = 0.379 (J0rgenson, Franx, & Kjaergaard 1992), although we find that its GC systems 
is better described with higher ellipticity, and use e = 0.6. For this galaxy, we do not need to model its GCs but can 
use the observed radial distribution of GCs from the WFPC2 observations of Harris et al. (2009). 

We then estimate the GC specific frequency, Sn, for each galaxy based on its absolute V magnitude, My. We use 
a linear interpolation of the Virgo early-type galaxy Sn data presented in Table 3 of Peng et al. (2008), giving us an 
estimate of the total number of GCs in each galaxy. Although we do not make any adjustments for morphological 
type, most cluster members in the core are early- types similar to the galaxies analyzed in Peng et al. (2008), and 
early-types also generally have higher Sn than late-types. 

Lastly, we determine an effective radius, R e , for the GC system, also assuming that it scales with galaxy My. We 
used a 2nd order polynomial fit between R e and My in which GC systems of more luminous galaxies are larger. We 
derived this relation using the GC systems of 100 Virgo cluster galaxies observed in the ACS Virgo Cluster Survey 
(Cote et al. 2004; Jordan et al. 2009). This survey observed early- type galaxies with a wide range of luminosities 
(—22 < Mb < —15). The spatial density profiles were fit with Sersic profiles to derive the effective radius of the GC 
systems as a function of galaxy luminosity. The details of this relation and its derivation will appear in a subsequent 
paper (Peng et al, in prep). The spatial extent of GC systems is generally larger than that of the stars, and is consistent 
with previous measurements showing that the local Sn for giant ellipticals increases with galactocentric radius. With 
this estimate of R e , we can then infer the contribution of GCs from galactic systems as a function of position in the 
central mosaic of observations. 

Using the same spatial grid shown in Figure |4j we estimate the total number of GCs expected in each cell from all 
of the galaxies in the Eisenhardt et al. (2007) catalog (not including NGC 4874). We correct this number for the total 
observed area — accounting for unobserved and masked area — and the local GCLF completeness. For each cell, we thus 
have an estimate of the total number of observed GCs that could be contributed by the GC systems of known cluster 
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galaxies. 

We then use this simulated distribution of GCs to statistically subtract the contribution from "galactic" GCs. In 
each cell, we sample a Poisson distribution about the expected value, and then assign each GC a random position 
in the cell. Additionally, we mask regions around all the known galaxies, as described above, with the results shown 
in Figure IA1I We do expect that some galactic GCs will contribute to the apparent intergalactic GC population, 
especially GCs from the other gE, NGC 4889, but the expected number and spatial distribution of the galactic GCs 
(i.e., concentrated around the galaxies) do not agree with the spatial uniformity and overall large numbers of GCs 
observed throughout the cluster core. 

We also note that although our modeling of GC systems only extends to galaxies with My < — 14 mag, this is 
unlikely to bias our final result. While it is true that Sn does tend to rise at lower luminosities, and some individual 
dEs have high Sn, the mean value is not very high. As shown by Lotz & Miller (2007) and Peng et al. (2008), dEs 
and dSphs have a very wide range of Sn, with some being high but many others quite ordinary or near zero. The 
mean Sn for dwarfs in the magnitude range that concerns us is ~ 2. Studies of Local Group dwarfs show that fainter 
than a certain luminosity (Mb ~ —10), dSphs do not appear to have GCs. 

A. 3. Systematic Errors From Modeling 

The technique we use to account for the unmasked galactic GCs introduces a systematic uncertainty into our analysis. 
We can quantify this effect by varying the assumed parameters over a wide range of acceptable values. We detail this 
in this section, but our main point is that our conclusions are not particularly sensitive to the specific parameters 
assumed for the modeling of galactic GC systems, because if Coma early-type galaxies have the same range of specific 
frequencies as Virgo early-type galaxies, then there are simply not enough galactic GCs to account for all the GCs 
that we observe. 

We have produced a range of plausible galactic GC populations, varying the Sersic n index of the assumed profiles, 
the effective radius, i? e , and the specific frequencies of the galaxies, Sn- By varying the Sersic index from n = 1.0 to 
n = 4.0 (exponential to de Vaucouleurs), the inferred number of IGCs varies from 47000 to 44000. The exponential 
profile is steepest in the outer regions and so gives the highest estimate for the IGC population. If we assume that all 
galactic GC systems have the shallower outer profiles of a de Vaucouleurs profile, then our inferred IGC population is 
smaller by ~ 3000 GCs. This effect is only at the peak-to-peak level of 6-7%, and is small enough to justify our choice 
of a single Sersic index for our GCS models. A single n — 2 model gives us an inferred population of 46500 IGCs, so 
we take the systematic error here to be (±3 ooo)- 

By varying the specific frequency, we affect the total number of galactic GCs. Our modeling uses the Sn — L 
relationship from Peng et al. (2008). Although there is no evidence that Sn for normal galaxies is different in the 
Coma cluster, we can make the assumption that Sn is systematically different from Virgo by ±20%. Doing so would 
change the number of IGCs by ±2, 000. This does not appreciably alter the significance of the detection. We include 
this in the systematic error budget. 

We can also test the dependence of inferred IGC numbers on the effective radii of the GC systems. We can assume 
that the R e of GC systems in Coma are a factor of two larger or smaller for a given galaxy luminosity than for their 
counterparts in Virgo (a fairly extreme assumption). If the GC systems are larger, the number of unmasked galactic 
GCs is increased, and the number of IGCs will go down. Using this assumption in our modeling, the total number of 
IGCs changes by ±3, 000. Once again, the overall effect on the significance of our detection is rather small. 

Given that our systematic errors are comparable to or larger than our Poisson errors (which are at the level of 
±1600 GCs), it is important to include these in the error budget. Nevertheless, it is clear that the details of the 
modeling does not affect the conclusion that there is a significant detection of IGCs in the Coma core. By combining 
the estimated errors independently and in quadrature, we estimate that there is a combined systematic error of 
±5000 IGCs - The inferred number of IGCs is thus 47000 ± 1600 (random) ± 4 5 °° a ° (systematic). 

A.4. A Test of GC System Modeling 

There are three outer ACS fields that we do not use for our background estimate because they are near larger Coma 
galaxies. Assuming that these fields have few or no IGCs, they present the opportunity to test the validity of our 
methodology. This test is not ideal, because NGC 4839 is a giant elliptical galaxy that is the dominant member of 
its subgroup, and our assumptions are most valid when we can average over many galaxies rather than be dominated 
by a luminous single galaxy. Nevertheless, we can apply the technique used to model the core region and check for 
consistency. Similar to NGC 4889, for which there are published measurements and no need to infer from scaling 
relations, we use the known Sn for NGC 4839 (Marin- Franch et al. 2002), although for every other parameter and for 
all other galaxies in the region, we use our scaling relations. We take the same approach to modeling and masking 
as was done in the core, although since the Eisenhardt et al. (2007) catalog does not cover this region, we use an 
equivalent cut with SDSS. For each field, after subtracting the modeled GC population and the estimated background 
from the total counts, we should expect to have numbers statistically equivalent to zero. For these three fields, we 
obtain residual counts of 13 ±9, 3± 11, and 10 ± 13, where the errors do not include the systematic errors as discussed 
above. We can assume that the systematic errors contribute roughly equally (or more) to the error budget. These 
numbers are satisfyingly consistent with zero, within the random and systematic errors. The predicted values do 
skew positive (with large errors) , but this is likely a reflection of our systematic error, and is consistent with what we 
discuss in Appendix IA.31 Over three fields, the residual is marginally significant over random errors, 26 ± 19 or when 
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corrected for the full GCLF, 0.007 ± 0.005 kpc~ 2 . This is still 8x smaller than our measured IGC surface density, so 
any systematic residual at this level would not affect the significance of our measurement. Even if we were to adjust 
our model to increase the predicted number of GCs so as to match these fields, the impact on our result and the IGCs 
numbers would be small and within our quoted errors. Thus, we are reassured that our GC subtraction procedure and 
estimation of our uncertainties is robust. 
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